Exact Asymptotics for One-Dimensional Diffusion with Mobile Traps 
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We consider a diffusing particle, with diffusion constant D' , moving in one dimension in an 
infinite sea of noninteracting mobile traps with diffusion constant D and density p. W e show that the 
asymptotic behavior of the survival probability, P{t), satisfies limt—,cx, [— In P{t)]/y^ p^Dt = A/^/n, 
independent of D' . The result comes from obtaining upper and lower bounds on P{t), and showing 
that they coincide asymptotically. We also obtain exact results for P{t) to first order in D' / D for 
an arbitrary finite number of traps. 

PACS numbers: 05.40.-a, 02.50.Ey, 82.20.-w 



The asymptotics of the survival probability, P{t), for 
a particle diffusing among mobile traps is a longstanding 
challenge. Over ten years ago Bramson and Lebowitz 
iQl showed by rigorous arguments that in one dimension 
P{t) has a stretched exponential decay for large t, 



P{t) - exp(-Ati/2) 



(1) 



The computation of the constant A, however, has thus 
far proved intractable, both analytically and numerically. 
This longstanding problem is resolved in this Letter. 

The problem addressed by Bramson and Lebowitz was 
first posed almost 20 years ago in the seminal paper of 
Toussaint and Wilczek (TW) These authors intro- 
duced the two-species annihilation process, A -f B 0, 
as a model of monopole-antimonopole annihilation in the 
early universe, though applications to chemical kinetics 
and condensed matter physics are more numerous 
TW showed that if the initial densities, /ca(0) and ps(0), 
of the A and B particles are equal (and the particles 
are randomly distributed in space), the densities decay 
asymptotically as i"''/^ in space dimensions d < 4. If the 
initial densities are different, however, the density of the 
minority species {A, say) decays much more rapidly. 

The connection with the trapping problem is as fol- 
lows. At late times pA{t) <C PB{t) and the A particles 
can be regarded as independently diffusing in a back- 
ground of the majority B particles, which act as traps 
for the A particles by virtue of the annihilation reaction. 
An equivalent problem, therefore, is to consider a single 
A particle moving among B particles (which do not in- 
teract with each other) and ask for the probability, P{t), 
that the A particle survives to time t. In the context 
of the original A + B ^ Q process {OT A + B ^ B [|, 
which has the same asymptotics [|l|) the "particle" A and 
"traps" B are taken to have the same diffusion constants, 
but for generality we will take them to be different, D' 
and D respectively. 

To our knowledge, no analytical result has been ob- 
tained up to now for the constant A in Eq. (^. Further- 
more, attempts to determine A by numerical simulations 
(or even to confirm the predicted stretched-exponential 
decay with exponent 1/2) are severely hampered by large. 



slowly decaying transients j^-^. In ref. 0, a sophisti- 
cated numerical approach enabled data to be obtained 
down to Pit) ^ lO^'^^, but still the asymptotic time de- 
pendence could not be unambiguously established. 

For this trapping problem, the traps are infinite in 
number, and distributed randomly on the interval (— oo, 
oo) with density p. By contrast, if the number of traps 
is finite, the problem is equivalent to the much studied 
predator-prey problem |p|,^ , where the traps are preda- 
tors and the particle is the prey. In this case the prey 
survival probability decays as a power law. 



P{t) i-S{Nl,Nr.D' /D) 



(2) 



where {Np) is the number of predators initially to 
the left (right) of the prey. To our knowledge the only 
exactly-solved examples are for Nl +Nji < 2. Obtaining 
analytical results for more than two predators is another 
longstanding challenge. 

In this Letter we obtain two exact analytical results. 
First, we finally resolve the question of the asymptotics 
of P{t) for the trapping problem. We verify the asymp- 
totic form (|l|), and determine exactly the value of the 
constant A, namely A — Ap{D /ttY/'^ , i.e. 



Pit) - exp[-4p(DV7r 



,1/21 



(3) 



Note that this asymptotic result depends only on the 
density, p, and diffusion constant, D, of the traps and 
is independent of the diffusion constant D' of the parti- 
cle. The value of D' , however, does affect the rate of 
approach to asymptopia. Eq. (||) is derived by obtaining 
upper and lower bounds for P{t), and showing that they 
converge for large t. Furthermore, the form of the lower 
bound suggests why large corrections to the asymptotic 
behavior might be expected. 

Secondly, we outline [ pd| a perturbation theory in 
D' /D, based on a path-integral formulation, for the com- 
putation of the exponent 9 in Eq. (|2|). For small D' / D 
and arbitrary values and Nn, we find 



O 



(4) 
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where N = + is the total number of predators 
and A = Nl — N^j measures their left-right asymmetry 
with respect to the prey. 

We will first describe the treatment of the trapping 
problem with infinitely many traps. The bounds that 
lead to Eq. (j^) are obtained as follows. 
1. Upper Bound. An obvious upper bound, Pu{t), on 
P{t), for any D' , is provided by the problem with D' — 0, 
in which the particle stays at its initial position, which 
we call a; = 0. Although we have as yet been unable to 
make this bound rigorous p^ , it is intuitively clear that 
when, as here, the traps are (statistically) symmetrically 
placed with respect to the particle, the particle will on 
average survive longer if it stays still than if it diffuses. 
This assertion has been checked, using the algorithm out- 
lined in ref. 0, for all (lattice) walks up to time t = 28. 
It is also supported by Eq. (^) , which shows that for any 
symmetric case (A = 0), the decay of P{t) is faster for 
small D' than for D' = 0, for any finite N. 

For D' — (sometimes called the "scavenger model" 
1^ in predator-prey terminology), P{t) is just the prob- 
ability that none of the moving traps has reached the 
origin up to time t. This problem is exactly soluble p4| . 
Since similar techniques will be needed to derive the lower 
bound, we outline the solution here. 

The A'^ traps move independently according to the 
Langevin equations 



(5) 



where r]i{t) is Gaussian white noise with mean zero and 
correlator 



(6) 



The quantity P{t) is just the product of the individual 
trap probabilities. For a given trap starting at Xi, the re- 
quired probability is Pi{xi,t) = eTf{\x^\/VlDt). So 
our upper bound is 



/ ^ 

Puit)^(Y[eTf(\x,\/ 



\i=l 




(7) 



where (. . .) means an average over the initial positions 
of the traps. Since the also independent, the lat- 

ter average also factors. Using N — pL, where L is the 
length of the system, and each Xi is uniformly distributed 
in (-L/2,L/2), gives 



Puit) 



1 - 1/L 



L/2 



L/2 

exp 



dxeTic{\x\/VADt) 



pL 



(8) 



where the final result follows on taking the limit L oo. 
2. Lower Bound. Consider the same system as before 
but with a pair of absorbing boundaries at x = ±Z/2. We 



consider the subset of initial conditions in which all the 
traps lie outside the interval [—1/2, 1/2) (and the particle 
is at X = 0), and trajectories in which neither the par- 
ticle nor any of the traps has crossed a boundary up to 
time t. We calculate the probability, PL^t), of such an 
occurrence over the ensemble of all initial conditions and 
trajectories. These restricted initial conditions and tra- 
jectories are a subset of all the possible initial conditions 
and trajectories in which the particle never meets a trap. 
It follows that P{t) > PL{t), i.e. PL{t) is a lower bound. 

The probability that there are no traps in the inter- 
val (—1/2,1/2) at t = is exp{—pl). Given that there 
are no traps in this interval at < = 0, the probability 
that no traps enter the interval up to time t is given by 
the same result as in the derivation of the lower bound, 
namely exp[— 4p(_Dt/7r)^/^]. Finally, the probability that 
the particle, starting at x = 0, has not left the interval 
{—1/2,1/2) up to time t is given, for times t ^ l^/D', 
by (4/7r) exp(— 7r^£>'t/P). Assembling these contri- 
butions gives 

P{t) > (4/7r) expi-ApiDt/TT)^/"^ - {pi + n^D't/l^)] . (9) 

Since this inequality holds for all I, the best lower 
bound is obtained by maximizing with respect to I. The 
optimum value is I = {27r'^ D't/p)^/"^, and the best lower 
bound is 

PL{t) = - exp[-4p(i?</7r)i/2 - 3(7r2p2i:)'t/4)i/3] . (iq) 

Since the second term in the exponent is negligible com- 
pared to the first as t — ^ cx3, the two bounds converge 
to yield the asymptotic form, Eq. (^), for P{t). More 
precisely, we can take the logarithm of P{t) and divide 
out the leading \/t dependence to get 



In P(t) 



< 



4 „ /TT 



TT^Vi jP'/D)^/^ 



(11) 



giving limt^oo -[liiP{t)]/{p^Dt)^/^ = 4/0F. 

As an aside we note that, while the left-hand inequality 
in ( p]) holds for all t, since Eq. (||) does, the right-hand 
inequality is strictly a large t result. This is because the 
factor (A/tt) exp{—Tr'^D't/P) in Eq. (||) comes from the 
lowest mode in the Fourier decomposition of the survival 
probability of the particle in the interval {—1/2, 1/2). This 
mode dominates for D't ^ P, which requires p^D't ^ 1. 
A lower bound on P{t) valid for all t can easily be writ- 
ten down by including all Fourier modes, but the laige-t 
form (^) is sufficient for present purposes. 

In Figure 1, the left and right hand sides of Eq. (12), 
representing the two asymptotic bounds, are plotted and 
compared to the numerical data of ref. Q . The data are 
plotted using the dimensionless time p^Dt, and the axes 
are chosen to test the asymptotic form (||). The data 
lie between the bounds (except at early times where lat- 
tice effects are important M). The bounds displayed are 



2 



for the case D' = D. The very slow convergence of the 
bounds is consistent with the observed trend in the data. 
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FIG. 1. Numerical data from ref. [7] and the upper and 
lower bounds, Eq. (12), for D' — D. The numerical data 
were generated with D — D' — 1/2 and p — 1/4. 



Some comments on Eq. (|11|) are in order. First, 
the right-hand inequality already suggests the possibil- 
ity that the approach to asymptopia could be very slow, 
as observed in simulations, since the second term van- 
ishes only as t~^^^. For the case D' — D that is usually 
simulated, the asymptotic result (^) gives, for a dimen- 
sionless time p^Dt = 10* say, a survival probability of 
order 10~^*. This is far smaller than can be reached 
in simulations (or experiments!) even with sophisticated 
methods 01 . Yet even at this large time, the second term 
on the right-hand side of (|ll|) is still 39% as big as the 
first! This term would give a further factor of 10"'^* in 
the survival probability, i.e. a total probability of order 
-[^Q-i36| Qj: course, this second term is only a bound, 
and the true correction to the asymptotic form could be 
smaller. The form of the bound is nonetheless suggestive 
and illustrates how large subdominant terms could arise. 
A clearer understanding of these subdominant terms is 
necessary before a detailed comparison with numerical or 
experimental data can be made. 

A second interesting point is that the mean-square dis- 
placement of the particle, averaged over the surviving 
trajectories used to derive the lower bound, grows as 
P ^ {D't/pY with V = 2/3. This can be compared 
with the estimate — 0.5 to 0.6 obtained from recent 
simulations 0. Note, however, that these simulations 
were not in the asymptotic regime, so the numerical esti- 
mates should be treated with caution. Furthermore, the 
value 1^ — 2/3 obtained from the lower bound for P{t) 
does not necessarily represent any kind of bound for v. 
Nevertheless, the rough agreement between the measured 
value and that obtained from our very simple arguments 



is again suggestive. 

The model can be generalized to include n diffusing 
particles starting from the same point M. Our bounds 
can easily be generalized to this case [Q. The proba- 
bility that all n particles survive at time t satisfies the 
same inequalities ( |TT| ) except that the final term on the 
right-hand side acquires an additional factor n^^^. So 
the bounds converge for t oo and the same asymp- 
totic form, Eq. (||), is obtained for all n. 

Our approach is also readily generalized to dimensions 
d > 1 111. For d < 2 we find P(t) ~ exp[-adpiDt)'^/% 
where = (2/7rd)(47r)^/2 sin(7rd/2) while, for d = 2, 
P{t) ^ exp[—4:Tr pDt / \nt]. The latter agrees with the 
functional form obtained in ||l[ , but with a precise value 
for the constant. For d > 2, simple exponential decay is 
obtained, in agreement with but the bounds no longer 
converge and the decay constant cannot be determined. 

In the remainder of this Letter we sketch ||ll[ the 
derivation of the perturbation theory, for a finite number 
of traps, that leads to Eq. (^). We recall that the N traps 
move independently according to Langevin equations (|^) 
with noise correlator (^ . We call the particle coordinate 
XQ{t). It obeys the Langevin equation Xq = 770 (t), where 
770 (t) is independent of the other noise terms, with corre- 
lator {r]o{t)r]o{t')) — 2D'S{t — t), corresponding to a par- 
ticle diffusion constant D' . It is convenient to introduce 
the relative coordinates yi — Xi — xq, i — 1, . . . , N. The 
corresponding Langevin equations are iji = S,i(t) where 
the noise terms, £,i{t) — rii{t) — 'rio{t), have correlators 



(12) 



where Dij = DSij 



D' . The equivalent Fokker-Planck 

equation, dP/dt — Y^f / ^Vi^Vji to be 
solved subject to absorbing boundary conditions at yi = 
for any i (i.e. when any trap meets the particle). In 
principle this equation can be solved by a coordinate 
transformation that diagonalizes the matrix D to give 
isotropic diffusion. This transformation, however, also 
rotates the edges of the absorbing region so that they are 
no longer mutually orthogonal but form the edges of an 
A^-dimensional hyperwedge. For N ~ 2 this becomes a 
two-dimensional wedge, and the problem can be solved 
exactly [^,0, but to our knowledge there are no exact 
solutions for TV > 3. 

An alternative starting point is the path-integral rep- 
resentation for the survival probability of the particle, 



Pit) = 



J^Dyjt) eMS[y\) 
jDyit) eM-S[y\) 



(13) 



where y 



IS a 



shorthand 

rt 



for (yi 



■■■iVn), S[y\ = 

i'^/'i)E"j=iiD'% lo dt' Ut')y3{t'), and the subscript 
R indicates that the path integral is restricted to "sur- 
viving" paths, in which none of the yi has changed sign 
up to time t. The matrix D is easily inverted to give 



3 



= (1/-D)(% - A), where A = D'/{D + ND') will 
be our expansion parameter. Thus 



S[y\^So[y\-SM 



1 - A 

-E, 



E fdt'Mt')? 
, "'0 



dt'Ut')yj(t') 



(14) 
(15) 

(16) 



where the diagonal terms, i — j, have been absorbed into 

Sq. 

A convenient normalization of the path integral is pro- 
vided by the path integral for the survival probability, 
Po{t) ~ t^^^^, of the same problem with D' — and 
D ^ Do = D/{1 - A). Then 



P{t) _ {exp{SM))i 



Po{t) (exp(5i[^)) 



(17) 



where the average in both cases is over paths weighted 
by exp(-S'o[2^). 

The final step is to expand numerator and denomina- 
tor using the cumulant expansion. To first order in A (i.e. 
first order in 5*1) one has 

P{t)/Poit)^exp{{Si)R~{Si)) . (18) 

Since the different values of i decouple in Sq we have 

{S^)n = {X/2D)J2 f dt' m'))R{yAt'))R ■ (19) 

The corresponding expression for the unrestricted paths 
vanishes, since {yi{t')) = by symmetry. The quan- 
tity {yi{t'))R is independent of i apart from the sign: 
{yi{i'))R > (< 0) for traps that stay to the right 
(left) of the particle. The calculation of {yi{t'))ji is 
straightforward. We only require the result in the regime 
to <S^t' t, where to = yfo/D and yio is the initial value 
of trap i. In this regime one finds (for traps which start 
on the right of the particle) ||ll[] 



{Ut'))R = 2iDo/Trt'y^^ , to « t' « t 



(20) 



We now insert this result into Eq. (|9|) . Let there be Nl 
(Nji) traps to the left (right) of the particle. Then the 
[NLiNL-l)/2+NR{NR-l)/2] pairs of traps whose 
members are on the same side of the particle contribute 
with positive weight to ( |l9| ) while the NlNr pairs whose 
members are on opposite sides enter with negative weight. 
Defining A = Nl - Nr (and N = Nl + Nr) this gives 



{Si)r 



XDq 
ttD 



(A'^-N) Int 



(21) 



to leading logarithmic accuracy for large t. The factor In t 
comes from the integral J^^{dt' /t') that appears when Eq. 



(gg) is substituted into Eq. (yj). Using the form (|2C|) for 
all t' , with t and to as upper and lower cut-offs, is cor- 
rect to leading logarithmic accuracy |l^. Note that the 
contributions from the lower cut-off, which depend on 
the initial positions of the traps, do not contribute to 
the leading logarithm. Using ( pi] ) in (p^, recalling that 
(Si) = and Po{t) ~ t^^^^, and noting that, to leading 
order in D' , A = D'/D and Dq = D, one finds P{t) ~ t'" 
with 9 given by Eq. (^ . The power-law decays hold when 
t > ylJD for all i. The magnitudes of the initial coordi- 
nates yoi determine the amplitude of the power-law, while 
their signs determine the exponent through the value of 
A. 

Eq. (3) can be compared with the exactly solved cases 
Nl= Nr = 1 ("surrounded prey") and A^l = 0, A^ii = 2 
("chased prey") §@, for which 9 = 7r/2a and [7r/2(7r- 
a)] respectively, where a = cos^^[D' /{D+D')]. Expand- 
ing these to first order in D' gives agreement with Eq. 
i). 

In summary, we have obtained the exact asymptotic 
form of the survival probability of a particle moving in 
an infinite background of mobile traps. The asymptotic 
form is independent of the diffusion constant of the parti- 
cle. A perturbation expansion in the ratio of the particle 
and trap diffusion constants has been developed for the 
case of a finite number of traps. 

We thank V. Mehra and P. Grassberger for supplying 
their raw data, and for a helpful correspondence. This 
work was supported by EPSRC grant GR/R53197. 
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